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Abstract In a recent work, Unruh showed that Hawking radiation is un- 
affected by a truncation of free field theory at the Planck scale. His analysis 
was performed numerically and based on a hydrodynamical model. In this 
work, by analytical methods, the mathematical and physical origin of Un- 
ruh's result is revealed. An alternative truncation scheme which may be more 
appropriate for black hole physics is proposed and analyzed. In both schemes 
the thermal Hawking radiation remains unaffected even though transplanck- 
ian energies no longer appear. The universality of this result is explained 
by working in momentum space. In that representation, in the presence of 
a horizon, the dAlembertian equation becomes a singular first order equa- 
tion. In addition, the boundary conditions corresponding to vacuum before 
the black hole formed are that the in-modes contain positive momenta only. 
Both properties remain valid when the spectrum is truncated and they suffice 
to obtain Hawking radiation. 
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1 Introduction 



The theory of black hole evaporation |l] is sore beset with two unsolved 
dilemmas: 

1 ) The transplanckian issue [[J : when deriving Hawking radiation in the 
usual framework of free field theory, one calls upon vacuum fluctuations 
at X~ |5] whose energies are 0(e M /M). These propagate as such up 
to a Planckian distance of the horizon where there energy is redshifted 
down to 0(1) and then further down to a typical frequency of 0{M~ 1 ) 
upon reaching 1 + . Since gravitational interactions become strong at 
the Planckian scale free field theory is at best dubious. 

2) The unitary issue : in the semi-classical theory of back reaction, both 
the matter which is the source of gravity (the star), and the "partners" 
to the emitted Hawking photons fall into the singularity, giving rise to 
a density matrix description of the radiation, thus, in the last analysis, 
to a non-unitary description of the evolution. Can one incorporate such 
a quandary into quantum mechanics? 

Both will require a deeper knowledge of how gravity reacts to Hawking emis- 
sion in order to be resolved. But perhaps at a more preliminary stage progress 
may be made by introducing an effective theory which one guesses incorpo- 
rates some of the features which might emerge from the fundamental theory. 
Such is the nature of a recent interesting contribution of W. Unruh 0, who 
addressed himself to the taming of the transplanckian monster. The present 
paper, inspired by Unruh's work, contains an analysis as well as a general- 
ization of the taming mechanism. 

Unruh's work is based on an analogy between the hydrodynamic equations 
of motion of sound waves in a moving fluid and those that governs s-wave 
emission of a massless scalar field from an incipient Schwarzschild black hole 
in free field theory. Thus one is led to predict the production of a thermal 
flux of phonons as the fluid passes from sonic to supersonic flow ||. Through 
numerical computation he subsequently showed |7|] that a truncation of the 
spectrum of sound -using a rather natural algorithm- in no way affected the 
thermal emission of Hawking phonons in the circumstances that this emission 
occurred in the untruncated theory. Carrying this lesson over to the black 
hole situation, the implication is that tinkering with the transplanckian part 
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of the photon spectrum will leave the thermal emission unaffected. We shall 
here follow up on Unruh's tinkering, first by introducing a truncation scheme 
which we believe to be more appropriate to the black hole situation, and then 
supply the mathematical rationale for the resistance of Hawking's result to 
such mutilation. Throughout we have preserved the linearity of the field 
equations. Whether non linear effects preserve Hawking radiation remains a 
moot point. 

The paper is organized in seven parts. Section |2| contains a brief re- 
view of Unruh's considerations. In preparation for our analysis we present in 
Section [3] a variant of the technique of Damour and Ruffini |9j , using momen- 
tum space considerations. This technique provides for a simple and elegant 
characterization of initial conditions which then rapidly leads to the under- 
standing which we seek. In Section |], the photon spectrum is truncated in 
this approach. It follows immediately from the formalism that tinkering with 
the transplanckian part of the spectrum does not affect the Hawking thermal 
emission. Section [5] incorporates the Damour- Ruffini technique into the Un- 
ruh hydrodynamic truncated model. Using similar reasoning as in Section 
|], the usual thermal emission is once more recovered. Section || contains a 
comparison of the wave packet trajectories in the two cases that have been 
analyzed so as to draw a physical picture of the production mechanisms. Fi- 
nally in Section [7| we speculate on the physics behind our phenomenological 
truncation procedure. Some interpretation in terms of quantum gravity is 
hazarded. The result of our analysis can be viewed as one of disappointment. 
Hawking radiation does not seems to be, in itself, an open door that leads 
to quantum gravity. Rather, it provokes thought in that direction, without 
offering, at least in a direct way, an orientation for the solutions. 

2 Unruh's Hydrodynamic Model 

The equation that governs the propagation of sound (amplitude = 0) in a 
perfect fluid in 1 + 1 dimensions which flows with a stationary background 
velocity field V(£) is 

[{d v + dsV(0) (d v + V(£)dt) - %]</> = . (1) 

This equation can be rewritten in d'Alembertian form □</> = with metric 
g 00 = 1, g 01 = g 10 = V, g 11 = — 1 + V 2 . It is readily diagonalized by the 
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transformation dt = drj + Vd£/(1 — V 2 ) ; d£ = dr to read 

aj ) =[(l-V 2 )- 1 d 2 -d r (l-V 2 )d r ](f) = . (2) 

In this form one recognizes the s-wave part of □ constructed from the 
Schwarzschild metric with the identifications: t,r = Schwarzschild time and 
radial coordinate respectively; V 2 = 2M/r. The horizon is at \V\ = 1. Thus 
a fluid, whose flow rate approaches V = — 1, will emit a thermal distribution 
of phonons whose temperature is given by (l/2n)dV/ dr\v=-i- At this point 
we simply use the form given by eq. ([!]) to motivate, with Unruh |7|], the 
truncation procedure used for this model. 

In the rest frame (V = 0), the spectrum of sound determined from eq. 
(HP is u — \k\. But since the fluid has atomic structure, the spectrum has the 
property to = \k\ when k < ko and ui — > ujq for k » ko- For example u = ko 
tanh (k/ko) might be expected to approximate the spectrum. In general for 
D = F(k), we are then led to modify eq. (JTp to 

[(-id v - idi-V) {-id v - Vid 6 ) - F 2 (-z^)]0 = (3) 

This equation has been the starting point of Unruh's numerical compu- 
tations. Propagating backwards in time an outgoing wave packet centered 
about a given negative frequency, he determines the Bogoljubov coefficients 
by decomposing the packet at early times into its positive and negative fre- 
quency components. Setting up a vacuum state at early times (See Section 
|B|), the P term in the Bogoljubov transformation (i.e. the weight of the posi- 
tive frequency part at early times of a negative frequency mode at late times) 
encodes the presence of outgoing phonons. This latter conforms to the ex- 
istence of the outgoing thermal Hawking flux - quite independently of the 
truncation function, F. 

An amusing aside is provided by Landau's theory of superfluid critical 
velocity, V s . Were vortices and rotons absent, the approach to critical super- 
fluid flow would be accompanied by a thermal shower of Hawking phonons. 

3 The Damour - Ruffini Method 

To explain and generalize Unruh's result we shall use the Damour- Ruffini 
technique 0, using a variant based on momentum space introduced in ref. 
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nj. It is especially convenient for our present purpose because the wave 
equation near the horizon and the boundary conditions defining in-modes 
take a particularly simple form. 

Equation (|2|) is valid outside a collapsing spherical star. First we trans- 
form it through use of the advanced Eddington-Finkelstein system, v and r, 
where v = t + r*,r* = r + 2M£n[(r - 2M) /2M] to give 

□0 = -[d r (l - 2M/r)d r + 2d v d r ]i/j(v, r) = (4) 

The solutions, 0, are connected smoothly to those which are inside the star. 
Since the details of the star's trajectory are irrelevant in the interesting 
asymptotic region where the star's surface approaches the horizon (r = 2M), 
we shall for simplicity of presentation take a very simple model: the star is 
defined by a thin shell radially falling inward with the speed of light. We 
choose its trajectory to be v — v st = where st labels the star's surface. 

Inside the star, for v < 0, the geometry is flat and described by the metric 
ds 2 = dv 2 — 2dvdr ( In terms of usual Minkowski time T one has v = T + r). 
In this coordinate system waves obey 

□0 = -d r (d r + 2d v )i) = (5) 
The usual Minkowski modes inside the star hereafter called in-modes, are 



xT {v \v) = e-^/V^X , (6) 
X r ( "V,r) = e-^-^/v 7 !^ . (7) 

The w-sector of the field operator is = / °° dX (a\xT^ + a \xT^*^) i- e - a 
superposition of in-modes given by eq. (0) inside the star. The Heisenberg 
state is the vacuum state |0j n ) which is annihilated by the destruction op- 
erators a\, i.e. there are no quanta inside the star. For a full description 
one should consider f-modes X\ as well. But these do not encode particle 
creation and will not be considered here. Henceforth we drop the label (u). 
(To be complete we also point out that in order to describe the s-wave sector 
of a 3 + 1 dimensional theory one should impose the additional boundary 
condition that the modes vanish at r = 0. This complicates the notations 
without modifying the result and will not be taken into account.) 
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Hawking radiation is encoded in the history of the outgoing modes xT as 
they propagate through the star into the space exterior to it. In particular we 
shall display their form outside the star both near and far from the horizon. 
In the first instance, the equation that governs their behavior is 

- d x [(x/2M)d x + 2d v ]iP{x,v) = (8) 

where x = r — 2M. This equation holds in the domain \xu\ << 1 where u 
is the eigenvalue of id v . Indeed the phase of the exact solution of eq. (^) 

e -iu>v e -i2u,(x+2Mlnx) &nd the phage q{ ^ so l ution Q f eq . (g) g-^g-^4Mlnx 

differ by 0(1) when \xu\ = 0(1). Henceforth we shall consider only typical 
uj = 0(M _1 ) which implies that x < M. 

From eq. (^), the w-modes (irrelevant for production) remain of the form 
e~ iuv / \/47ru , whereas the w-modes are given by the linear combination 

T^v) = ^==[6{x)Ax AMiuJ + 9{-x)B\x\ 4Mw } (9) 

Standard Klein Gordon normalization prescribes \A\ 2 — \B\ 2 = 1. 
Far from the horizon (x ^> 2M) , the wave equation reduces to 



-d a! {d x + 2d v )ijj = (10) 

and the -u-modes ipZ ut associated to asymptotic quanta on X + are equal to 
^out _ e -iu{v-2r) i y / 4^J_ Since the exact solutions of eq. (^) which connects 
to this asymptotic form is ip° ut = e - luJ ( v ~ 2r *) j *J the out-modes near the 
horizon are given by A = 1, B = in eq. (0): 

= -^=e~ iu)V 6{x)x iMiu} . (11) 
V Anu 

As the crux of our analysis lies in a careful formulation of the matching 
conditions at the star's surface which define the modes outside in terms of 
the modes inside the star we now go into this matter in some detail. In order 
to determine the out-particles content of the in-vacuum, it is propitious to 
compute the Fourier transform of the modes at the surface of the star (v = 0). 
The matching conditions are then implemented in simple and elegant fashion. 

The Fourier transform of eq. (§) is 

[pd p + 4Mioo + 1}^(P) = (12) 
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whence 



= [Cp-^9(p) + D\p\- iMi »9(-p)\ . (13) 

The domain of validity of eqs. [T2]) and is \p\ ^> M _1 (for typical u = 
0{M- 1 )). 

The Fourier transform of the Minkowski modes eq. (|7|) are proportional 
to S(p — 2A), and we recall that positive norm modes have A > 0. 

The object of the exercise is to use continuity to change the basis from 
modes inside the star (= xTip)) to those outside the star which for large 
p are of the from J-^(p). Since the former set have p > (since A > 0), 
continuity prescribes that p > is valid for the latter set as well. In this way 
one establishes that an equivalent set of in modes is found from Eq. ([13]) with 
D = and C = 1 (this set was obtained independently from Damour-Ruffini, 



but by working in Kruskal coordinates by Unruh |l Ij and Hawking [g] ) . We 
call this bases i>™(p) 



C» = ^(P)?^ 1 • (14) 

This gives ?p™(p) for sufficiently large p (e.g. p ^> M" 1 ). Thus for these 
values of p one finds that the expansion coefficients a^x which give ip l ™{p) in 
terms of x\ n (p) are (M / 8ti 2 A) 1//2 (2X)~ iMluJ . For smaller values of p which are 
relevant for the construction of packets which cross the surface of the star 
at large values of x, as we have stated above, ip™{p) and xT(p) coincide. In 
this case a^x — > S(lj — A). 

To see explicitly that only the large values of p are relevant for packets 
which reach X + at late times, one builds a wave packet with the modes eq. 

(ED: 

r -(uj-ui ) 2 /2a 2 

= ]du> (27r)1/4ff i/ 2 e-^CCp) (15) 

The phase e tuJU ° centers the wave packet on the light ray v — 2r* = uq and the 
gaussian e"^"^ l 2a centers the frequency around uq. Inserting the form of 
ip™ one obtains 



rJn (.. _\ _ V r M(27r) 1/4 (T 1 / 2 1 ( t , + 4Mlnp-Tto) a <7 3 /2 r iu (v+4M Inp-Ug) f-,a\ 
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whereupon it is seen that the wave packet is centered on the trajectory v + 
4M\np = uq. The momentum p at the surface of the star v = is given by 
p = e"°/ 4M . Hence for wave packets which reach X + at late times Uq, p at 
the surface of the star is much larger than M _1 . 

In conclusion, we have the important result that, for all values of p, 
only positive ones appear in the basis functions ip™(p). Furthermore, for 
p ^> M^ 1 which are those necessary to describe outgoing packets which 
begin their journey from the surface of the star to X + at values of x s t <C M, 
the content of the X l \(p) modes in terms of ip l ™(p) contain both signs of uo. On 
the contrary for p — > 2u, which are the relevant values to describe outgoing 
packets which begin their journey to X + at values of x s t 3> M, uo becomes 
equal to A. Therefore it is the former set that gives rise to particle creation 
on X + , whereas the latter give rise to no creation. In this formalism, this 
is what expresses the well known fact that Hawking radiation sets in as the 
star's surface approaches the horizon exponentially closely. To derive these 
results all that has been used is continuity at the star's surface. From now on 
we shall characterize in- vacuum by a^lOm) = where the annihilation 

operators associated with the ip™ modes and we re-emphasize that u takes 
on both signs in this characterization. 

To find the content of |0j n ) in terms of the out-modes (defined on X + ) one 
considers the Fourier transform of ip in (p) valid for p ^> M" 1 , hence x -C M 



V2tt Jo 



Y(-AMiuj)\x\ AMiul [e(x)e + ^ Mu> + 6(-x)e- 2nMw ] . (17) 



V2- 
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The mode ip™ lies on both sides of the horizon. Only the piece outside the 
horizon (i.e. proportional to 9(+x)) propagates out to X + . From eqs. (|TT|) 
and (0) one has outside the horizon {x > 0) 



tf(«)«(-H.) = { Sgt) " <°0 < 18 ' 
where we have introduced the Bogoljubov coefficients 
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O' 



7T 

Y(+AMiu)^2Mu e ~ 27TMuJ 



/7T 



l/SLKI = e" 4 ^ . (19) 

This implies a steady thermal flux of emitted particles at temperature T# = 
1/8ttM. 



The condition p > used in eq. (14]) in the present formalism is equivalent 



to the Damour-Ruffmi requirement that ip™{x) be analytic in the upper half 
x plane. What we have seen here is that it is a direct consequence of the 
existence of vacuum in the star. The description we have given here is readily 
adapted to take into account truncation of the transplanckian spectrum (p > 
1). 



4 Confronting the Transplanckien Hiatus 

The mechanism of how Hawking photons get created is through the combined 
gravitational and Doppler red shift which is incurred as a wave packet voyages 
from small values of x to I + . In the steady state of radiation the star's surface 
is exponentially close to the horizon (x st — Me~ l ^ 2M ). In consequence the p 
values, obtained from eq. fllTD , which dominate the integrand in the packet 
in this region are exponentially large: the saddle point of the integrand in 
eq. (0) is at p* = AMlu/x Jig. 

To see this more precisely, recall that the locus of saddle points traces 
out the classical trajectories of photons. Near x = they are obtained from 
the Hamiltonian constraint 



H = p (H + 2pv ) (= 0) (20) 

where we have used Eddington-Finkelstein coordinates x,v and their conju- 
gate momentum p and p v (compare with eq. (§)). The canonical equation 
are 

p = —dH/dx ; x = dH/dp ; v = dH/duj ; p v = , (21) 
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where the dot denotes derivative with respect to an affine parameter along 
the trajectory. These combine to yield 

p v = —uj = const, 
dp p 
dv ~ 4M 

On mass shell, H vanishes and accordingly 

4Mu; 

x = 

p 

Hence 

p(v)=p 8t e- v '* M and x(v) = ^^ e ^ 4M (25) 

Pst 

where p st is a constant of integration. From this last equation, we see that 
a photon of fixed energy u, erected at late time and reaching x = O(M) at 
t>* (where t>* is typically of order the life time of the black hole = 0(M 3 )) 
must have crossed the surface of the star at x st = 

with an enormous momentum p st = e +v *^ M 0{uj). This is the transplanckian 
problem. At these high values of p, and concommitantly small values of 
x, free field theory certainly breaks down. The gravitational interaction of 
the s-wave modes under consideration, both with other modes and with the 
background field whose source is in the degrees of freedom of the star becomes 
enormous. Therefore the underpinnings of the free field calculation become 
completely fallacious. 

This does not mean, however, that the derivation given in Section |3| is 
completely lost. This is Unruh's main point. We shall first illustrate these 
considerations in the context of the Damour-Ruffini formalism of Section |3] by 
introducing a truncation in that scheme and in Section |5] proceed to analyze 
Unruh's effective theory as given in Section 0. In both of these schemes one 
explores the hypothesis that the physics near the horizon can be mimicked 
by modifying the way the matter field propagates, but leaving both gravity 
and the linear character of the fluctuations unaffected. 

First a few words of clarification are in order. Any truncation scheme can 
be formulated in intrinsic geometric terms. However it is convenient to work 
in a coordinate system which is privileged in the geometry of the incipient 
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(24) 



black hole. We make the assumption that the truncation takes a simple form 
in such a privileged system. Explicitly, the origin is fixed and the angular 
momentum expansion is carried out with respect to it. The radius has an 
intrinsic geometric meaning. Furthermore the space time outside the star is 
static and the origin of time is given by the inception of radial infall. 

Let us therefore begin by truncating in the Eddington-Finkelstein system. 
We suppose that the energy spectrum gets modified for modes probing space 
time at scales less then 1, i.e. for p > 1. Then eq. fl8|) is assumed to be 
modified to the form 



g{-id x 



x 



2M 



g{-id x ) - 2id v 



i[)(x, v ) = 



(26) 



Hence the mode equation corresponding to Eq O) becomes 



[d p g(p) + 4Miu]r u (p)=0 . (27) 

The physics that goes into the evaluation of g(p) will be the subject of the 
conjectural discussion of Section [j]. For the present we only require that 
g(p) = p for p < 0(1) and g{p) = 1 for p > 1 since we anticipate that it 
is only the Planckian and transplanckian modes which will be affected by 
gravity. As in Unruh's model, the function g(p) is now to be identified with 
the energy associated with the mode number p rather than p itself. Indeed 
in flat space the proposed modification eq. (p6|) reads 



g(-id x M-id x )-2id v ]ip = (28) 

so that one sees that for outgoing modes A, the eigenvalue of id v , is equal to 
g{p)/2 rather than p/2 and is bounded by 1/2. In the curved space outside 
the star, the energy at X + (= uj) is Doppler and gravitationally red shifted 
when compared to the energy measured by the free falling observer (for a 
discussion of the red shift see ref [0 Chapter 3 and Appendix D). This red 
shift is represented by the factor x/4M in eq. (|26|) and the energy measured 
by the free falling observer is AMuj/x = g(p) which is now bounded by 1/2. 

It is instructive to see how the truncation deforms the classical trajectories 
from the geodesic associated with the free field. The truncated version is 
described by the Hamiltonian 
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n 



9(p)[w779(p) + 2Pv] (= 0) 



(29) 



and eqs (|^) and (|24D now read 



Pv 

dp/dv 



x 



-g(p)/4M 
4MLu/g(p) 



—uj 



(30) 
(31) 
(32) 



Integrating, one sees that \p\ decreases with v, initially (when \p\ > 1) linearly 
and then exponentially (when \p\ < 1). For \p\ > 1, one has \x\ = \4Mu\ 
and for \p\ < 1, \x\ is proportional to |p| _1 , hence increases exponentially 
with v. Thus at the same time as one avoids transplanckian energies, one 
ceases to approach the horizon to within transplanckian distances (typically 
ujM = 0(1)). Section ^| contains a sketch of the production process based 
on these classical trajectories (see Fig. |j). 

To establish Hawking radiation in the truncated theory we once more 
have to characterize the in-modes exterior to the star. Subsequently it must 
be shown that these in-modes evolve towards X + so as to give the required 
radiation. 

From the truncated equation inside the star (eq. (pg|) ) it is seen that 
positive A implies positive p. Thus continuity across v = once more implies 
that in-modes on the outside also have positive p. (This crucial result, which 
is independent of the details of g(p), results from our assumption of postu- 
lating a linear field equation which keeps the d'Alembertian factorized into a 
u and v piece. That this is sufficient is obvious, but that it is not necessary 
will become apparent in Section ^). Thus our in-modes, solution of eq. (p7|) 
with p > are 



The Fourier transform of ip l ™(x) has contributions from p > 1 for x = 0(1) 
and p < 1 for x > 0(1). This is seen by examining the saddle point condition 
in the Fourier transform which in fact reconstructs the classical trajectories 




(33) 
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eq. (|32|). The piece which is relevant for Hawking production is at x » 1, 
hence concerned with cisplanckian values of p and one has 



thereby recovering eq. QT7DH. 

It should be noted that the norm of the cisplanckian modes is the same as 
in the free field case. This results from the fact that the truncated problem 
is possessed of a conserved norm (= 2ti Jdp ip^ip) g{p) i>u>'(p) — <K<^ ~ ^0) 
which reduces to the free field form for \p\ < 1. Thus the decomposition of 
V^ n ( x )lcis ^ n ^° Pi eces localized at x > and x < as in eq. ([H]) has its usual 
interpretation in terms of unitarity. 

When translated into the variable x, what we have shown is that the 
free field vacuum is maintained at x > 0(1) since one only requires the 
characterization of the in-state in terms of its cisplanckian content in this 
region. Most succinctly, Unruh vacuum is still a valid concept just outside a 
Planckian skin of the horizon. This is every bit as valid as the statement that 
Unruh vacuum is the correct description of the in-state in free field theory 
outside the star. A little bit of Planckian fuzz around the horizon does no 
injury to the physics since the conversion to Hawking photons on I + occurs 
outside the star due to the redshift which the outgoing modes feel on the 
scale of x = O(M) and not x = 0(1) (this has been pointed out by Jacobson 
[ |j~3|j who however did not derive the insensitivity of Hawking radiation to 
planckian tinkering). 

1 We point out that eq. ( p8| ) has in general a large number of linearly independent 
solutions (an "infinity" if g(d x ) is not polynomial in d x ). These linearly independent 
solutions are recovered from the function ip™(p) by specifying which integration contour 
is used to take the Fourier transform. According to how the singularities of ip™(p) are 
encircled, different solutions are obtained. We have chosen the contour to coincide with 
the path in the absence of truncation, namely to lie on the real axis. Other prescriptions 
would presumably give rise to runaway solutions and the properties of the theory at low 
momentum would not coincide with the free field theory. 




(34) 
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The picture that emerges is that fluctuations within this skin steadily 
develop into outgoing pairs. Note that from eq. ([H]), p in this region grows 
linearly in v. The interpretation is that there is a conversion of modes within 
the skin to become the usual outgoing modes of free field theory, thereby 
guaranteeing a steady state. Thus at large values, p has become converted 
from an energy to a mode counting parameter. Note that the total time of 
evaporation is 0(M 3 ), so that the total number of modes which boil off from 
those initially inside the skin is 0(M 2 ) which is proportional to the usual 
estimate for the entropy of the black hole. 



5 The Truncated Unruh Model 

In this section we analyze the production of Hawking phonons based on 
Unruh's truncation eq. @. As in Section f|, the analysis is based on the 
momentum representation of modes. Complications occur because the equa- 
tion for the modes near the horizon is second order and it is necessary to 
control that the transplanckian sector does not contaminate the cisplanckian 
physics. This was obvious in Section ^ once the boundary conditions were 
set. 

To have a clear idea of the mechanism of production from the modes, it 
is first worthwhile to go into the classical motion along trajectories. Whilst 
this part of the analysis does not give the production per se, it does give the 
motion of wave packets before and after production. In this way one has a 
guide to the portion of the mode analysis which is relevant. 

The Hamiltonian, which generates the classical trajectories corresponding 
to the wave equation eq. (0) is 

H= [p v +pV(Z)] 2 -F 2 (p)] (=0) (35) 

In the context of eq. (|3~5"D, p v is the momentum conjugate to rj; p being 
conjugate to £. From the canonical equations one finds 



Pr! = —uj = const. , (36) 

dp/dn = - P V'(0 , (37) 

if- ii j V(0 + F'(p) on trajectories of type II or III , , 

1 V(£) — F'(p) on trajectories of type I 
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Since the fluid flows to the left with V decreasing to the left (so that near 
the horizon V = — 1 + £/4M), it is seen that \p\ always decreases in 77. As a 
consequence we can, and shall, describe evolution with respect to 77 in terms 
of p. 




Figure 1: The two curves —uj ± |.F(p)| (solid curves) and the family of lines 
—V(C,)p (doted lines) are plotted as functions of p for some representative 
values of £. Their intersections give the trajectories £(p). Their are three 
trajectories labeled I, II, III. The set I corresponds to a v = const trajectory. 
The set II crosses the horizon when £ = and then starts to propagate for 
£ < O(-AMuj) whereupon it corresponds to a u = const trajectory. The 
points in class III never reach the horizon: there are no solutions in this 
class for £ < £ m i„(u;) = 0(AMuj) but their are two solutions for £ < £ m j n 
corresponding to a trajectory which bounces. 

To analyze the trajectories we use the on mass shell condition — pV(£) = 
—uj ± In Fig. ^] is plotted as functions of p, the family of curves 

—V(£)p, and the two curves — uj ± The intersections give the trajec- 

tories at fixed u. It is seen that there are three types of trajectories: 

I The sequence of points near the origin all have small negative p, hence 
correspond to a usual null geodesic in Schwarzschild geometry which, 
in E F coordinates is v— constant. They are therefore uninteresting for 
the production of particles. 
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II These are points which, for £ ^> — 4Mu, have p < 1, hence F(p) = 
1. A particle on this part of the trajectory is at rest with respect to 
the fluid. For £ < —AMu, p is < 1 and the particles then propagate 
with respect to the fluid. This trajectory stops at p = owing to the 
monotonous character of p[rj) ( eq. (j37j)). 

Ill The third trajectory has p > 0. Once more for large positive p it 
is non propagating. As \p\ decreases, it starts to propagate for p = 
0(1) whereupon it reaches a minimum value of £(= £ m m(^)) which is 
0(AMu). As |p| continues to decrease, £ now increases so as to describe 
the propagating part of the trajectory for negative values of p. 

Trajectories II and III are plotted in Fig|4| where the functions £(t>) are 
displayed (these are very similar to the functions £(77) since the coordinate 
transformation used to go from eq. (P to eq. (|2|), shows that near the 
horizon 77 differs from v by a regular function of £.). Production is concerned 
with the mixing of trajectories say III into II, i.e. a wave packet localized 
on trajectory III at large £ and large p (that is at early times, rj small) has 
an amplitude (3 to be localized on trajectory II at small p (rj large). Note 
that had we taken to < rather than uj > 0, trajectory II would then have 
positive p and trajectory III negative p rather than the reverse. 

We now study eq. (|J) near the horizon where V ~ — 1 + £/4M. In 
momentum space, after some elementary manipulation, the equation takes 
the form 



p 2 



iuAM + 1 . a ' t 

d p + + 1AM 

p 



AM 2 

For \p\ < 1 where F(p) = p, the solutions are 



+ F 2 (p) U w (p)=0 . (39) 




I : UP) = \l—0(±p)\ P \- 4muJ - 1 e- l8Mp ; \p\ < 1 



II, III : UP) = \-0{±p)\p\~^- 1 ; \p\<\ . (40) 

7T 

Solutions in class I are associated with the trajectories of class I (v-modes) 
and are not involved in production (the singular behavior of these modes is 
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due to the presence of the (spurious) Cauchy horizon at £ = 8M) . For p 1 , 
where F(p) = 1, the solutions are 



^(p) = ^9(±p)\p\- 4M ^- 1 e +l4Mp \p\ ia± ; p»l (41) 

where a± are the roots of a(a — 1) — (4M) 2 = 0. We take M ^> 1 so that 
a± ~ ±4M — |. These modes serve as a basis of second quantization. Their 
positive and negative frequency parts will then correspond to the bases of the 
quantized field wherein the annihilation part (positive frequency) annihilates 
the in- vacuum (i.e. vacuum at large £, hence large p). To determine this we 
refer to the conserved scalar product which we norm in the conventional way 



i J d£ [<M0* (d v + v(Od^) MO - MO* <- MO] = ±<J(w - u') 

(42) 

or in momentum space 

i2tx J dp [M(pT (d v -ip- 1/4M - p/AMd p ) M?) 

-M(pY ^ 4>M\ = ±6(u — u r ) . 

(43) 

Inserting eq. ( ^T|) into eq. ( ^3|) it is seen that the sign of the scalar product 



is determined by the sign of a± in eq. (0). The condition of positive 
frequency (a+ = +4M) corresponds to taking wave packets localized along 
the trajectories — pV{^) = —oj + |F(p)| rather than — pV(£) = —uj — \F{p)\. 
The additional requirement that the modes be localized along trajectories II 
or III rather than I then imposes p > 0. The condition of positive p which 
played an essential role in the preceding sections is thereby recovered. In 
summary the in-modes appropriate for second quantization which give rise 
to Hawking radiation are 

4%(p) = J— gr p } p -*Miu-l e -i±M Pp i4M+l/2 . p>1 (44) 
V 7T 

It remains to establish how these in-modes which have been characterized 



at large p evolve into their forms eq. (4"0D at small p. To this end we use 
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the WKB approximation whose validity we justify subsequently. Thus we 
approximate the full solution of the in-mode for positive p by 



iiM\p-J P dp F(p)/p] 

ti?(p) = e{p) P -' m - 1 r _= — . (45) 

lF{p)/p 



Equation ^g) is exact for both \p\ < 1 and \p\ ^> 1. The validity of the WKB 

Jity 

< 1 . (46) 



expansion is assured owing to the inequality 

d(p/4MF(p)) 



dp 

This condition can be given a geometric interpretation. To this end one 
reexpresses the solution £ (p) of the on mass shell condition H = ( eq. (j35|) 
as C(p) = £f(p) + A£(p) where £f(p) corresponds to a point which moves 
with the fluid to + pV(£p) = and A£ describes motion with respect to 
the background (pV'(£ F )A£) 2 - F 2 (p) = 0. Then the validity of the WKB 
approximation takes the form 

— -— ^ — — < 1 . 47 

dp 

which expresses that the motion with respect to the background fluid is suffi- 
ciently "slow". The l.h.s. of this inequality is 0(M~ r ) which by assumption 
is 1. Any back scattering due to the second order character of the differ- 
ential equation will be a non perturbative effect (typically of 0(e~ M )) which 
would result in the mixing of u and v type modes. It is negligible. 

Once having established that the low p part of the in-modes is of the form 
eq. ( |40|) , types II and III, with only positive p, the results of Section |] follow 
forthwith: to wit Fourier transform gives a thermal distribution in u arising 
from packets built from the low momentum part of the modes. 

In any truncation scheme, be it linear or nonlinear, one may hope that 
a condition of adiabaticity similar to eq. (ff7|) will be applicable. It will 
then imply that the creation of particles at a Planckian distance from the 
horizon is strongly suppressed. Thus Unruh vacuum will once more be a valid 
concept on scales x > 0(1) and the usual spectrum of Hawking radiation will 
be recovered. Were there any Planckian particles created at x = 0(1) they 
would severely disrupt the Hawking flux. Their absence can be taken to 
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be an expression of the general principle (very well verified experimentally) 
that flat space is stable against creation of Planckian particles. Since the 
curvature in Schwarzschild geometry acts on scales Ax = O(M), at a Planck 
distance from the horizon space looks almost flat and the principle can be 
applied with confidence. Thus the absence of created Planckian particles at 
the horizon should come as no surprise. 



6 The classical trajectories 




Figure 2: The classical trajectories of outgoing null geodesies (thin dot- 
ted curves) given by the Hamiltonian eq. are displayed in the Edding- 
ton-Finkelstein coordinate system. The light ray which generates the horizon, 
the infalling shell, the origin r = and the singularity are also represented. 

We now have on hand three schemes on how to get Hawking radiation. 
It is interesting to display in an Eddington-Finkelstein graph the classical 
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trajectories corresponding to the modes in the three cases. These encode 
vacuum fluctuations in the past which are converted into quanta at x = 
0{M). 

In Fig H is displayed the usual free field model. A vacuum fluctuation 
emerges from the star, it is a pair that straddles the horizon. The fluctuation 
hugs the horizon at an exponentially small distance. Outside the star it starts 
to propagate (outwards for the observed Hawking photon and towards the 
singularity for its partner -unobservable in Schwarzschild time but taking 
form on the other side of the horizon in finite Kruskal time). Upon reaching 
x = O(M) the fluctuation has become an on mass shell quanta which now 
propagates along the trajectory v — 2r = const. 

In Fig |3] the effect of the simple truncation of Section ^ is shown. We 
have taken g(p) to be unity for p > 1 (and of course g(p) = p for p < 1). The 
propagation on either side of the horizon begins at the edges of the Planckian 
zone |x| = 0(1). The Planckian zone |x| < 1 is thus steadily solicited to give 
out radiation at a steady rate. 

Finally in Fig. |] the trajectories in Unruh's truncation are sketched. 
Here they begin from I~, but unlike the free field case they do not go into 
the star to reflect and then come out as in Fig |2|. Rather they reach the 
horizon region directly whereupon at x = ±1 they reflect. The amplitude 
of the reflected wave is augmented; it is accompanied by a partner which 
appears on the other side with amplitude (3 and the production is encoded 
in a 2 — 1 = f3 2 = {er* — It is fairly difficult to interpret this structure 

in terms of the black hole. In the fluid it comes about due to the "riding in" 
of the modes on the background flow. 

More important, from the figures it is seen that in all cases, at large dis- 
tance from the horizon \x\ 3> 1, \p\ <C 1, one recovers the usual free particle 
trajectories. Therefore, independently of the truncation we have used, Hawk- 
ing radiation remains a pair production phenomenon: the outgoing quanta 
are accompagned by partners on the other side of the horizon. 

7 Conclusion 

It appears to us that the scenario based on eq. (|26| ) could well turn out to 
reflect some part of the truth. The function g(p) deviates from its free field 
value, p, for p > 1 so as to give rise to a sort of quasi-particle description. 
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particle 



partner 




Figure 3: The trajectories of outgoing light rays given by the truncated Hamil- 
tonian eq. \2$ ) are displayed in the Eddington-Finkelstein coordinate system. 
In this case the trajectories no longer approach the horizon exponentially but 
rather they stick at a planckian distance. 

The strong gravitational interaction among the modes will result in contin- 
uous mixing of angular momenta, so in the rigorous theory, restriction to 
s-waves will no longer be possible. The whole medium is to be regard as 
a matter gravitational soup, which has some s-wave content. For values of 
x greater than unity this content becomes that of the usual free field and 
there should be a turn over from a quasi to true particle theory. The trans- 
planckian soup steadily feeds into the free field sector. It should be possible 
to display the transition region (x = 0(1), p = 0(1)) by perturbative meth- 
ods, for example by expanding the gravitational action to quadratic terms 
in fluctuations around the background geometry. This will show how the 
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modes start interacting with each other as they move into the Planckian re- 
gion. We might expect g(p) to become complex for p > 1, corresponding to a 
Planckian lifetime of the s-wave quasi-particle. Such behavior might encode 
that fact that an s-wave gets swallowed up in the extrapolation backwards in 
time towards the Planckian skin. Then instead of the extrapolation through 
the star drawn in Fig |3], the modes just peter out within the skin. This is 
indicated by some shading in Fig [5]. 

In addition there is the interaction of the modes with the gravitational 
field emanating from the sources which constitute the star (possibly taken 
together with the mean effect due to the other modes, as in the semiclassical 
theory). Here it would be "recoil" effects of the gravitational field which 
would be responsible for deviations from free field theory. That such effects 
can be important has recently been shown by one of us (RP) in the context 
of the accelerated detector [TJJ], and more recently corroborated in a study 



of accelerating mirrors ||15|| . What the incidence of these two types of effects 
on the modes is for the unitarity problem remains to be seen. 

Finally one must be prepared to encounter the very strong coupling prob- 
lem which will arise well into the transplanckian region where conceptual 
problems will arise concerning the nature of space-time. 

Whatever, on the basis of the above considerations, we conjecture strongly 
that Hawking radiation is protected from the vicissitudes of quantum gravity. 
It appears as an essentially kinematic response to the presence of the event 
horizon encountered in the collapse of a macroscopic black hole. 



Added note After this manuscript was completed, C. Bouchiat and F. En- 
glert called our attention to the following conceptual problem. The trun- 
cation we have used treats u-modes and v-modes asymmetrically, thereby 
explicitly breaking the invariance of the theory under local Lorentz transfor- 
mations. However one should recall that the formation of a black hole by 
the collapse of a star induces an asymmetry between u-ness and v-ness: only 
a future horizon is formed but no past horizon exists. Concommitantly the 
field state is asymmetric, the v-part being Schwarzschild (Boulware vacuum) 
in character whereas the u-part is Kruskal (Hartle-Hawking) in character (i.e. 
the state is Unruh vacuum). 

We beleive that it is not unreasonable to envisage that a phenomenological 
truncation of the field equations is state dependent since it should encode 
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the dynamics of the matter field state induced by the quantum gravitational 
interactions. In Unruh vacuum the truncation would then treat u and v 
modes differently. What would be the effect of a symmetric truncation in a 
symmetric state such as the eternal black hole situation remains to be seen. 
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particle 




Figure 4: The trajectories of light rays in Unruh's truncated hydrodynamic 
model are displayed in the Eddington-Finkelstein coordinate system. Only 
class II and III light rays have been displayed as it is these which are respon- 
sible for production. Note that since v and 77 differ by a regular function of 
x, the trajectories in the coordinate system rj, £ given by the Hamiltonian eq. 
\3dj ) would be very similar to those depicted in the figure. 
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Figure 5: The trajectories of light rays in the Eddington-Finkelstein coordi- 
nate system if the truncation g(p) where complex. In this case the trajectories 
disapear into some quantum fuzz which is represented by some shading in the 
figure. 
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